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ABSTRACT 

Relations between the observed quantities for a beamed radio jet, apparent transverse speed and apparent 
luminosity (/3 app ,L), and the intrinsic quantities, Lorentz factor and intrinsic luminosity (j,L ), are investigated. 
The inversion from measured to intrinsic values is not unique, but approximate limits to 7 and L a can be found 
using probability arguments. Roughly half the sources in a flux density-limited, beamed sample have a value 
of 7 close to the measured (3. dpp . The methods are applied to observations of 119 AGN jets made with the 
VLBA at 15 GHz during 1994-2002. The results strongly support the common relativistic beam model for an 
extragalactic radio jet. The (/3 app ,L) data are closely bounded by a theoretical envelope, an aspect curve for 
7 = 32, L a = 10 25 WHz" 1 . This gives limits to the maximum values of 7 and L in the sample: 7 max k, 32, 
andL„. raax - lO^WHz" 1 . No sources with both high /3 app and low L are observed. This is not the result of 
selection effects due to the observing limits, which are flux density S > 0.5 Jy, and angular velocity /1 < 4 mas 
yr" 1 . Many of the fastest quasars have a pattern Lorentz factor 7^ close to that of the beam, 7$, but some of 
the slow quasars must have j p <C Three of the 10 galaxies in the sample have a superluminal feature, with 
speeds up to /3 app rj 6. The others are at most mildly relativistic. The galaxies are not off-axis versions of the 
powerful quasars, but Cygnus A might be an exception. 

Subject headings: BL Lacertae objects: general — galaxies: active — galaxies: individual: Cygnus A — 
galaxies: jets — galaxies: statistics — quasars: general 



1. INTRODUCTION 

In recent years VLBI observations have provided many ac- 
curate values of the apparent luminosity, L, of compact radio 
jets, and the apparent transverse speed, /3 app , of features (com- 
ponents) moving along the jets. These quantities are of con- 
siderable interest, but the intrinsic physical parameters, the 
Lorentz factor, 7, and the intrinsic luminosity, L„, are more 
fundamental. In this paper we first consider the "inversion 
problem;" i.e., the estimation of intrinsic quantities from ob- 
served quantities. We then apply the results to data from a 
large multi-epoch survey we have carried out with the VLBA 
at 15 GHz. 

The inversion problem is discussed in ^2]-^4]with an ideal- 
ized relativistic beam, one that has the same vector velocity 
everywhere, and contains a component moving with the beam 
velocity. The jet emission is Doppler boosted, and Monte- 
Carlo simulations are used to estimate the probabilities asso- 
ciated with selecting a source: that of selecting (/3 app ,L) from a 
given (7,L ); and the converse, the probability of (j,L ) being 
the intrinsic parameters for an observed (/3 app ,L). 
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In $4]we introduce the concept of an "aspect" curve, defined 
as the track of a source on the (/3 app ,L) plane (the observation 
plane) as 6 (angle to the line-of-sight, LOS) is varied; and an 
"origin" curve, defined as the set of values on the (7X0) plane 
(the intrinsic plane) from which the observed source can be 
expressed. These provide a ready way to understand the in- 
version problem, and illustrate the lack of a unique inversion 
for a particular source. Probabilistic limits provide constraints 
on the intrinsic parameters for an individual source; but when 
the entire sample of sources is considered, more general com- 
ments can be made, as in ^7] 

The observational data are discussed in §|5] They are 
from a 2-cm VLBA survey, and have been published 
in a se ri es of papers: IKellermann et al.l (1 19981 hereafter 
iPaperfl). IZensus et al.l (120021 [Paper IJ). IK ellerma nn et al.1 
d2004llPaper Ilj) iKovalev et al.l d2005llPaper IVh , and E. Ros 
et al. 2007, in preparation. This is a continuing survey and the 
speeds are regularly updated using new data; in this paper we 
include results up to 2006 September 15. The analysis also 
includes some results from the MOJAVE program, which is 
an exte nsion of the 2-cm surve y using a statisticall y complete 
sample dLister & Hom an 2005). Prior to lPaoer III}, the largest 
compi lation of interna l motions was in IVermeulen & Cohenl 
d 1994 hereafter lVC94l) . who tabulated the internal proper mo- 
tion, n, for 66 AGN, and /3 app for all but the two without 
a redshift. The data came from many observers, using vari- 
ous wavelengths and different VLBI arrays, and consequently 
were inhomogeneous. The 2-cm data used here were ob- 
tained with the VLBA over the period 1994-2 002, and c om- 
prise "Excellent" or "Good" apparent speeds ( Paper nj) for 
components in 1 19 sources. This is a substantial improvement 
over earlier data sets, and allows us to make statistical studies 
which previously ha ve not been possible . Other recent sur- 
veys are reported by iJorstad et al.l (l2005l hereaf ter |J05|) with 
data on 15 AGN at 43 GHz. lHoman et al.l (1200 lb w ith data on 
12 AGN at 15 and 22 GHz. and lPiner et all (l2006h with data 
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on 77 AGN at 8 GHz. 

In some sources it is clear that the beam and pattern speeds 
are different, and to discuss this we differentiate between the 
Lorentz factor of the beam, 7b, and that for the pattern, 7 p . In 
most of this paper, however, we assume 7b ~ 7 P and drop the 
subscripts. In $6]and §|7]peak values for the distributions of 7 
and L„ in the sample are discussed, and in ^8]the low-velocity 
quasars and BL Lacs are discussed. It is likely that some of 
these have components whose pattern speed is significantly 
less than the beam speed. The radio galaxies in our sample 
are discussed in $9] and we show that most of them are not 
high-angle versions of the powerful quasars. Cygnus A may 
be an exception, and we speculate that it contains a fast central 
jet (a spine), with a slow outer sheath. 

In this paper we use a cosmology with Hq = 
70 kms -1 Mpc -1 , ft m = 0.3, and fl\ = 0.7. 

2. RELATIVISTIC BEAMS 

In this Sectio n the standard relations f or an ideal relativis- 
tic beam (e.g., [Blandford & Konig|11979j) are reviewed. The 
beam is characterized by its Lorentz factor, 7, intrinsic lumi- 
nosity, L„, and angle 9 to the line of sight. From these the 
Doppler factor, 5, the apparent transverse speed, /3 app , and the 
apparent luminosity, L, can be calculated: 

( 5 = 7 - I (l- / 3cos6>r I , (1) 
~ ^ (2) 

L = L S n , (3) 

where f3 = (1 -7~ 2 ) I//2 is the speed of the beam in the AGN 
frame (units of c) and L is the luminosity that would be mea- 
sured by an observer in the frame of the radiating material. 
The exponent n in equation [3] combines effects due to the K 
correction and those due to Doppler boosting: n = a+ p, where 
a is the spectral index (S y v a ) and p is the Doppler boost ex- 
ponent, discussed in 32.11 

From equations (Q]i and ©, any two of the four parameters 
/3 app , 7, S, and 9 can be used to find the others; a convenient 
relation is /3 app = f]jS sin 9. Figure QJ shows S and /3 app as 
functions of sin 9, all normalized by 7; the curves are valid 
for 7 2 > 1. When sin 6* = 7" 1 , (5 = 7 and /3 app = /3 app , max = 
/3j. The "critical" angle 9 C is defined by sin# c = 7"', and the 
approximation 9/9 c ks 7 sin (9 will be used; this is accurate for 
7 2 > 1 and 9 2 < 1, and is correct to 20% for 9 < 60° and 
f3> 0.5. 

It is also useful to regard /3 app and 9 as the independent 
quantities. Figure [TJ) shows S and 7 as functions of sin#, all 
normalized by /3 app . The curves are calculated for /3 app = 15 
and change slowly with /3 app , for /3 app ^> 1. 

2.1. Assumptions 

In our analysis we assume that a source contains an ideal 
relativistic beam; one that is straight and narrow, and in which 
the pattern speed is the same as the beam speed: 7 P ss 7b. In 
particular, the Doppler factor for the core must derive from 
the same values of 7 and 9 as apply to the value of /? app for the 
moving component, several pc or more away. Many sources, 
however, are seen to have more than one moving component, 
and they may have different values of /3 app . In these cases we 
have selected the fastest speed, on the grounds that of all the 
components, it is the one most likely to be moving at near 
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7sin9 /S app sin0 

FIG. 1. — Top: Parameters for a relativistic beam having Lorentz factor 7 
and angle to the LOS 8. (a) Curves are plotted for 7=15 but change slowly 
with 7, provided 7 2 2> 1. (b) Curves are plotted for /3 app = 15 and change 
slowly with /3 app , provided /3^p p S> 1. In (a) the quantities are normalized by 
the constant Lorentz factor; in (b), by the constant apparent speed. Bottom: 
Results from a Monte-Carlo simulation of a flux density-limited survey se- 
lected from the parent population described in Appendix [A] (c) Probability 
density p(6|7f) and cumulative probability P(6> |7f) (heavy line) for 7 fts 15. 
Roughly 75% of the selected sources will have 7 sin 6* < 1; i.e., 8 < 8 C . Val- 
ues of 7 sin (9 < 0.15 and > 2.0 are unlikely; the cumulative probabilities are 
approximately 0.04 and 0.96, respectively, (d) Probability density and cumu- 
lative probability for selecting a source at angle 8, for /3 app ~ 15. As /3 app 
decreases the probability curve becomes more peaked, and the peak moves to 
the left. 

the beam speed. It probably is due to a shock associated with 
an outburst in flux density, wh ile some of the slow er compo- 
nents might be trailing shocks (lAgudo et alJ feOOl). The main 
shock itself must be moving faster than the beam, but the syn- 
chrotron source, which is a density concentration behind the 
shock, can have a net spe ed slower than the shock, as shown 
by numerical simulations ( Agudo et al. 200j]). 

In some sources the only component we see is stationary 
at a bend in the jet. We believe that in these cases we see a 
standing shock, or perhaps enhanced radiation from a section 
of the jet which is tangent to the LOS. These components will 
have 7 P <C 7b and might be part of the population of slow 
quasars discussed in $8] 

It is clear that some jets are not straight, and that 9 is not 
the same in t he core and in the moving components. See, 
e.g., 3C 279 dHoman et al.l l2003h where the velocity vector 
changed during the course of observations, and 0735+178 
and 225 1+158, where the image shows a jet with sharp bends 
(Paper I). However, in cases where at least moderate superlu- 
minal motion is found, the motion must be close to the LOS, 
and any changes in angle will be strongly amplified by pro- 
jection. An observed right angle bend could correspond to an 
intrinsic bend of only a few degrees. 

The Doppler boost exponent p depends on geometry an d 
optical depth and is discussed by iLind & Blandford! d!985l) . 
For a smooth jet, p = 2, and this value is appropriate for the 
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core region, where, we assume, the relativistic beam streams 
through a stationary t = 1 region. In some cases the mov- 
ing component can be modelled as an isolated optically thin 
source, which would have p = 3, but for most of the sources 
the flux density is dominated by radiation f rom the core, and 
we use p = 2 here. See Figure 5 in lPaper IVl The other term in 
the exponent n is the spectral index a. Nearly all the sources 
have a "flat" spectrum, |a| < 0.5, and they also have variable 
flux density and a variable spectrum. Because of the time de- 
pendence it is not possible to generate a useful index for each 
source, and we take a = as a rough global average, giving 
n = 2. This is further justified in 35.31 where it is shown that 
n = 3 does not fit the data. However, this choice of a = 
clearly leads to error for those sources with a high Doppler 
factor, say 5 = 30, because the K-correction must cover a fre- 
quency range of a factor of 30. This introduces uncertainty 
into the estimates of intrinsic luminosity. 

We shall use equation (0 as if L is independent of 9, but 
this is not necessarily so. The opacity in the surrounding ma- 
terial may change with 9, and the luminosity of any optically 
thick component may change with angle. Other changes in L„ 
might be caused by a change in location of the emission re- 
gion, as 9, and the refore the Dopp ler factor and the emission 
frequency, change (Lobanov 1998). 

3. PROBABILITY 

The probability of selecting a source with a particular value 
of 9, 7, /3 app , or S from a flux density-limited sample of 
relativistically-boo sted sources is central to our discussion. 
Because S oc S 2 ( §5.31 and S decreases with increasing 9, 
the sources found will preferentially be at smal l angle s, even 
though there is not much solid angle there. IVC94I calcu- 
lated the probability p(9\ ^f) (the subscript f mean s fixed) in 
a Eucl idean universe, and lLister & Marscherl d 19971 hereafter 
LM97) extended this with Monte-Carlo calculations, to in- 
clude evolution. However, the observations directly give /3 app , 
not 7, and p(9\fi ap p.f) is generally not an analytic function. 
To deal with this, M. Lister et al. (in preparation) are using 
Monte-Carlo methods to study the probability functions. We 
use one of their simulations here, as an illustrative example. 

In the Monte-Carlo calculation, a simulated parent popu- 
lation is created (see Appendix [A), from which one hundred 
thousand sources with S > 1 .5 Jy are drawn. We select a slice 
of this sample with 14.5 < 7 < 15.5 and form the histograms 
in Figure [TJ;, showing the probability density p{9\~fs) and the 
cumulative probability P(9\ji) for those sources with 7 w 15. 
The histograms vary slowly with 7, provided j 2 ^> 1. They 
are similar to the equiva lent diagrams calculated by VC94 
(Figure 7) and by LM97 (Figure 5). Figure QJ may be di- 
rectly compared with Figure QJ, which is a purely geometric 
result from equations Q] and [2] The peak of the probability 
is at sin# « O.6/7, where /3 app k, O.97 and 5 ks 1.5j. The 
50% point of P(9\jf) is at 7 sin 6* w 0.7, giving a median value 
0med « O.77- 1 « 0.79 c . 

An interesting measure of the cumulative probability is 
P(9 = 9 C ), the fraction of the sample lying inside the critical 
angle. The slow variation of this fraction with 7 is seen in 
Figure [2k a rough value is 0.75; i.e., most beamed sources 
will be inside their "I/7 cones." In this paper we will take 
0.04 < P < 0.96 as a practical range for the probability. This 
corresponds, approximately, to 0.15 < 9/9 c < 2 for 7 = 15, 
and the angular range for this probability range varies slowly 
with 7. Figure l9l (Appendix lAl shows the (7,6*) distribution 
for 14,000 sources from the simulation, along with the 4% 
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FIG. 2. — (a) Probability that 7sin$ < 1 as a function of 7. (b) (solid 
curve) median value of 7/ /3 app and (dashed curve) median value of /3 app sin 8, 
as functions of /3 app . 

and 96% limits. 

We have now described p{9\~fi), the probability for selecting 
a jet at angle 9 if it has Lorentz factor 7/. However, given 
that we observe /3 app and not 7, we must consider also the 
probability ;?(f?|/3 app .f); i.e., the probability of finding a jet at 
the angle 9 if it has a fixed /3 app . We again use a slice of the 
Monte-Carlo simulation, now for 14.5 < /3 app < 15.5, to get 
the probabilities shown in Figure [TJl. The probability density 
curve is broad, and as /3 app decreases it becomes more peaked. 
The median value of /3 app sin 9 is shown with the dashed line 
in Figure^, as a function of /3 app . 

The probability p(7|/3 a pp,f) is a ls° of interest. Figure [3^ 
shows an example, for /3 app rs 15. The probability is sharply 
peaked at 7 ~ /3 app . The median value is 7med/Aipp = 1 -08, and 
it changes with /3 app as shown with the solid line in Figure^. 
The sharp peak can be understood in geometric terms. In Fig- 
ure QJ) one sees that there is a large range of 9 over which 
7 changes little from its minimum value near /3 app , and Fig- 
ureQJl shows that most of the probability is in this range. For 
about half the sources with /3 a , 



15, 7 is between 15 and 
16, but the other half is distributed to 7 = 32, as shown in Fig- 
ure [3^. For lack of better information, it often is assumed in 
the literature that 7 w /3 app , but this is not always valid. 

Figure [3p shows p(S\f3. dpvS ) and P(<5|/3 apPjf ) for /3 app « 15. 
The curves change slowly for /3 app 3> 1. Unlike the Lorentz 
factor, the probability for the Doppler factor does not have a 
sharp peak. Consequently L a , which varies as S 2 , is poorly 
constrained by /3 app . 

In this paper a particular Monte-Carlo simulation is used, 
to show probability curves in Figures Q] and [3] and, numeri- 
cally, to find the 4%, 50% and 96% levels of the cumulative 
probability distributions. These are fairly robust with regard 
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FIG. 3. — Probability density and cumulative probability (heavy line) when 
P m « 15. (a) p( 7 |/3app,f); (b) pi5\P mf ). 



to evolution and parent luminosity functions. We have com- 
pared them among several of the simulations calculated by 
M. Lister et al., and the variations are not enough to materi- 
ally affect any of the conclusions in this paper. 

Figures QJ[3] are not valid in the non-relativistic case, where 
1 < 1, 7 « 1, S « 1, Aw w /3sin0, and p(9) ~ sin<9. Our 
discussion is also not valid for samples selected on the basis 
of non-beamed emission. 

4. THE INVERSION PROBLEM 

VLBA observations can directly give apparent speed /3 app 
and apparent luminosity L, but the Lorentz factor 7 and the 
intrinsic luminosity L are more useful. We refer to the esti- 
mation of the latter from the former as the inversion problem. 

The inversion is illustrated with Figure |4] On the left is the 
intrinsic plane, with axes 7 and L„, and on the right is the 
observation plane, with axes /3 app and L. Consider a source at 
point a in Figure gji, with 7 = 20 and L„ = 2 x 10 24 WHz" 1 . 
Let it be observed at 9 = 1 .3°, so that /3 app = 15.0 and L = 2.2 x 
10 27 WHz" 1 . This is the point z in Figure |4j?. Now let 9 vary, 
and the observables for source a will follow curve A. We call 
A an aspect curve. It shows all possible observable (/3 app ,L) 
pairs for the given source a. The aspect curve is parametric in 
9, with 9 = on the right, as shown. The height of the curve 
is fixed by the value of 7, and the location on the x-axis is 
fixed by 7 and L . The width of the peak is controlled by the 
exponent n in equation[3] as discussed in ^5.31 

Now consider a source with observational parameters at 
point z in Figure What can be said about the intrinsic 
parameters for this source? From equations QJ[3] curve Z in 
Figure^ can be drawn; Z contains all possible pairs of intrin- 
sic parameters from which source z can be expressed. We call 
Z an origin curve. It is parametric in 9, with 9 = on the left 
as shown. The curve has been truncated at 7 = 32, because 
this is the approximate upper limit of 7 for our data, as shown 
in m 

Given the lack of a constraint on 9, the inversion for the 
observed point z in Figure|4]3 is not unique. Any point on the 
origin curve Z in Figure^ could be its counterpart. This gives 
limits to 7 and L„, but they usually are broad. The limits get 
tighter when the probability of observing a boosted source is 
considered, as in the next section. More general results apply 
in a statistical sense when a sample of sources is considered. 
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FIG. 4. — Illustrating the intrinsic (left) and the observation (right) planes 
for relativistic beams. The origin point a in panel (a), with 7 = 20, can be 
observed anywhere on the aspect curve A in panel (b), by varying 8. The 
observed point z in panel (b), with /3 ap p = 15, can be expressed from any 
point on the origin curve Z in panel (a). Both curves are parametric in 8 with 
8 increasing as shown. The maximum of the aspect curve in panel (b) is at 
8 = 2.9° and /3 a pp = 19.97. The minimum of the origin curve in panel (a) is 
at 8 = 3.8° and 7 = 15.03. Panel (c): as in panel (a) but with the origin curve 
truncated at points g and h, the 4% and 96% cumulative probability limits, 
respectively. Panel (d): as in panel (b) but with the aspect curve truncated at 
points u and v, the 4% and 96% cumulative probability limits, respectively. 



4. 1 . Probability Cutoffs 

The probabilities associated with observing beamed 
sources were discussed in §|3] We now use the 4% and 96% 
cumulative probability levels to define the regions where most 
of the sources will lie. Figures 0J; and[4]i repeat Figures [4^ 
and|4]3 with the origin curve truncated at P(7|/3 app ,f) = 4% and 
96%, and the aspect curve similarly truncated at P(9\jf) = 4% 
and 96%. Note that points g and h do not correspond to points 
u and v. The probabilities can be seen in Figures [TJl and QJ, 
respectively, as functions of sin 9. 

The luminosities are double-valued in Figure |4] The prob- 
ability cutoffs are found by integrating along curves A and Z, 
and not by accumulating values of 7 or /3 app along both sides 
of the minimum, or peak. An example of accumulating 7 on 
both sides of the minimum of an origin curve is in Figure [3^. 

In Figure |4j; the points on Z have different Lorentz factors, 
but all have /3 app =15. The run of 7 vs 9 along curve Z is 
shown in Figure [5] which essentially is a section of the curve 
in Figure [TJ). The probability p{9\(3 s , vv s), shown in Figure [TJl, 
varies slowly along this curve, and is indicated with the line 
width. 

From Figure 2J; we now have probabilistic limits for the 
intrinsic parameters of the observed source z. Points g, h, 
and the minimum give 15 < 7 < 25.6 and 1.0 x 10 24 W 
Hz" 1 < L a < 5.4 x 10 25 W Hz" 1 . Note that these values do 
not describe a closed box on the (7,L a ) plane. Rather, the 
possible values must lie on curve Z. The highest 7 goes with 
the lowest intrinsic luminosity, the lowest 7 goes with an in- 
termediate luminosity, and the highest luminosity goes with 
an intermediate 7. 
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1 2 3 4 5 6 
Aspect Angle (cleg) 

FIG. 5. — Curve Z from Figure[4j; is shown on the (7, 9) plane. The prob- 
ability density for rinding a source with /3 a pp ~ 15 is indicated by the width 
of the line. Points g and h are the locations where the cumulative probability 
f (0|Aipp.f) reaches 4% and 96%, respectively. For a range of 9 around the 
maximum probability at 9 ~ 2.5°, the value of 7 changes slowly. As shown 
more directly in Figure[3k approximately half the sources with /3 ap p = 15 will 
have 7 between 15 and 16. However, 9 is not similarly constrained. 

A large survey will likely contain other sources with /3 app 
near that of source z. They will have various luminosities 
and will form a horizontal band in Figure 0Ji. That group 
of sources will have a distribution of 7 with minimum 7 m j n w 
/3 app and median 7 me d ks 1.1 /3 app , according to Figure^. This 
means that, for any individual source, it is reasonable to guess 
that 7 is a little larger than /3 app , although that guess will be far 
off for some of the objects. It is correct to say that about half 
the survey sources with /3 app s=s 15 will have 15 < 7 < 16, and 
that about 95% of them will have 15 < 7 < 25.6. The value 
95% results from the 4% above point g in FigureHj;, and 1% 
above 7 = 25.6 when the curve is continued above point h. 

5. THE DATA 

The 2-cm VLBA survey consisted of repeated observations 
of 225 compact radio sources, over the period 1994-2002 . 
Since that time the MOJAVE program dLister & Homanl2005T) 
has continued observing a smaller but statistically complete 
sample of AGN. Most of the sources have a "core-jet" struc- 
ture, with a compact flat-spectrum core at one end of a jet, 
and with less-compact features moving outward, along the 
jet. The VLBA images were used to find the centroids of the 
core and the components, at each epoch, and a least-squares 
linear fit was made to the locations of the centroids, relative 
to the core. The apparent transverse velocity was calcul ated 
from the angular velocity and the redshift. See iPaper IIll and 
E. Ros et al. 2007, in preparation, for details. 

Each component speed is assigned a quality factor Excel- 
lent, Goo d, Fair, or Poor according to criteria presented in 
IPaper lilt but only the 127 sources with E or G components 
are used here. Eight of the 127 are conservatively classified 
by us as Gigahertz-Peaked-Spectrum (GPS) sources. This 
classification is given only to so urces that have always met 
the GP S spectral criteria given bv lde Vries. Barthel. & O'Deal 
(119971) . and is based on RATAN mo nitoring of broad-ba nd in- 
stantaneous radio spectra of AGN (Kovalev et al. 1999) 10 . In 
GPS sources the bulk of the radiation is not highly beamed, 
as it must be if our model is to be applicable, and we omit the 
GPS sources from this study. The final sample contains 119 
sources, comprising 10 g alaxies, 17 BL Lac ob j ects, and 92 
quasars, as classified by I Veron -Cettv & Veron (2003). (See 
classification discussion in IPaper IVh The sample and the 

10 See also spectra shown on our w eb site 
http : / /www .physics .pur due . edu/astro/MOJAVE/ 



(/3 app ,L) values used here are given on our web site'O. The 
/3 app data are updated from values in Paper III with the addi- 
tion of results from more recent epochs given in E. Ros et al. 
2007, in preparation, and on the web site. 

Values of (/3 app ,L) for the 119 sources are plotted in Fig- 
ure [6] Error bars are derived from the least-squares fitting 
routine for the angular velocity. The luminosities are calcu- 
lated, for each source, from the median value of t he "total" 
VLBA flux densities, over all epochs, as defined in Paper IV. 
SvLBA,med is the integrated flux density seen by the VLBA, 
or the fringe visibility amplitude on the shortest VLBA base- 
lines. The luminosity calculation assumes isotropic radia- 
tion. Error bars are not shown for the luminosities. Actual 
erro rs in the m easurement of flux density are no more than 
5% (IPaper IV)) . but most of the sources are variable over time 
(see lPaper IVl Figure 11). 

An aspect curve for 7 = 32, L = 10 25 WHz" 1 is shown in 
Figure [6] It forms a close envelope to the data points for L > 
10 26 W Hz" 1 . At lower luminosity the curve is well above the 
data, and, as shown in $7] and $8] lower aspect curves should 
be used there to form an envelope. A plot similar to the one 
in Figure [6] is in IVermeulenl (1995), for the early data from 
the Caltech-Jodrell Bank 6-cm survey (lTavloretal.lfT 996). 
Although no aspect curve is shown in IVermeulenl d 1 9951) . it is 
clear that the general shape of the distribution is similar at 6 
and 2 cm. The parameters of the aspect curve in Figure [6] are 
used in $7]to derive limits to the distributions of 7 and L a for 
the quasars. 

5.1. Selection effects 

A striking feature of Figure [6] is the lack of sources to the 
left of the aspect curve; i.e., we found no high-/3 app , low- 
L sources. We recognize two possible selection effects which 
might influence this, the lower flux density limit to the survey, 
and the maximum angular velocity we can detect. We now 
combine these to derive a limit curve. 

The 2-cm survey includes sources stronger than 1 .5 Jy fo r 
northern sources, and 2.0 Jy for southern sources (IPaper J) . 
Additional sources which did not meet these criteria, but were 
of special interest, are also included in the full sample. How- 
ever, here we a re using the median VLBA flux density values 
from IPaper IVl for the sub-sample of 119 sources for which 
we have good quality kinematic data, and the median of these 
values is 1.3 Jy. We choose S m \ n = 0.5 Jy as the lower level of 
"detectability," although 10% of the sources are below this 
limit. The completeness level actually is higher, probably 
close to 1 .5 Jy, but the survey sources form a representative 
sample of the population of sources with SvLBA,med > 0.5 Jy. 

The angular velocity limit, /i max , is set by a number of 
factors, including the complexity and rate of change of the 
brightness distribution, the fading rate of the moving com- 
ponents, and the interval between observing sessions. These 
vary widely among the sources, and there is no easily quan- 
tified value for ^t max . In practice, we adjusted the observ- 
ing intervals for each source according to these factors, with 
AT being about one year in most cases. This was usually 
sufficient to eliminate any ambiguity in defining the angular 
velocity as seen on the "s peed plot s," the position vs. time 
plots shown in Figure 1 of Pa per IlJ . For some sources there 
was little or no change in one year, and these were then ob- 
served less frequently. For others, a one year separation was 
clearly too long, and they were observed more frequently, typ- 
ically twice per year for complex sources. The fastest angular 
speeds we measured were < 2 milliarcsec (mas) yr" 1 , and we 
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FIG. 6. — Values of apparent transverse speed, /3 a pp, and apparent luminosity, L, are plotted for the fastest E or G component in 119 sources in the 2-cm 
VLB A survey. The aspect curve is the locus of (/3 a pp , L) for sources with 7 = 32 and L a = 1 X 10 25 WHz 1 , as 8 varies. Curve K is an observational limit set at 
SvLBA.med = 0.5 Jy and fx = 4masyr~' ; the hatched region is usually inaccessible. The horizontal lines are the minimum values of redshift, ZminO'^appX for which 
the angular velocity is belo w the limit, < 4 mas yr~' . The vertical lines are the maximum values of redshift, Zmax(L), for which the flux density is above the 
limit, S > 0.5 Jy. See j|5.1| Red open circles are quasars; blue full circles, BL Lacs; green triangles, galaxies. 



saw no evidence for faster motions that would require more 
frequent observations. It is important to note that even pro- 
grams with shorter sampling intervals, down to every 1 or 2 
months, have not detected many speeds o ver 1 mas yr" 1 , and 
none significantly larger than 2 mas yr " 1 dGomez et al.ll200il 
iHoman et al.ll2001l: Uorstad et alj|2005l) . 

A rough limit on our ability to identify very fast compo- 
nents is given by our typical one year observing interval and 
the fading b ehavior of jet comp onents. From an analysis of 
six sources, IHoman et al.l (120021) found that the flux density 
of jet components fades with distance from the core as R~ 13 . 
If a jet component is first identified at a separation of 0.5 mas 
with a flux density of 50 mJy, that component will probably 
have faded from view when it is 4 or 5 mas away, where it will 
have a flux density of only a few mJy. Such a component, ap- 
pearing just after a set of observations, could fade from view 
before the next observation a year later, if it was moving at > 4 
mas yr" 1 . In practice, however, we would be likely to observe 
such a source in the middle of its cycle, and it would appear to 
have jet components a few mas from the core which flicker on 
and off in an unpredictable fashion. So while we would not 



have been able to measure the actual speed of such a source, 
it would have been identified in our sample as unusual, and 
followed-up with more frequent observations. Given that we 
identified no such objects, we take 4 mas yr" 1 as a reasonable 
upper limit to the speeds we are sensitive to with our program. 

It is possible that some components could fade more rapidly 
than the above estimate, and if so, our limit would have to 
be reduced accordingly. There is some evidence that rapid 
fading occurs at 43 GHz, and in $8] we describe a source with 
a component moving more rapidly at 43 GHz than at 15 GHz. 
It is likely that the difference is due to a combination of a fast 
fading rate and better angular resolution, combined with the 
shorter observing intervals, at 43 GHz. Even here, however, 
the observed speed at 43 GHz is well under the limit. 

Curve K in Figure [6] parametric in redshift, is calculated 
from the limits S = 0.5 Jy and fj, = 4.0 mas yr" 1 . The hatched 
region to the left of the curve is inaccessible to our obser- 
vations except in special circumstances, such as when the 
brightness distribution is simple and there is only one feature 
in the jet. The horizontal lines in Figure [6] show the mini- 
mum redshift associated with a value of /3 app , set by the dis- 
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FIG. 7. — As in Figure |6]but with curves for 3 values of n, the Doppler boost exponent. The curves all have 7 = 32, and are truncated at the 4% and 96% 
cumulative probability limits. Values of L are adjusted to optimize the fit near the top and the right-hand side. 



tance at which /1 = 4 mas yr ; while the vertical lines show 
the maximum redshift associated with a value of luminosity, 
set by the distance beyond which the flux density is below 
0.5 Jy. Thus, every point to the right of curve K has a range 
of redshift within which it is observable, and that range fixes 
a spatial volume. Inspection of the diagram shows that the 
volume goes to zero at the limit curve and increases towards 
the envelope. This gradient constitutes the selection effect. 
Sources are unlikely to be found near the limit curve because 
the available volume is small. The volume inreases towards 
the envelope; and, for example, at L = 10 26 W Hz" 1 , /J ap p = 20, 
the range 0.08 < z < 0.30 is available. In the sample of 1 19 
sources we use, there are 10 sources in this range, all of them, 
evidently, far from the region in question. At L = 1 25 W Hz" 1 , 
/3 app = 10, the range 0.04 < z < 0.10 is available and 6 of the 
survey sources are in this range; again, none of them is near 
the region in question. Hence, the lack of observed sources 
to the left of the envelope is not a selection effect; but rather, 
must be intrinsic to the objects themselves. 

5.2. The Fast Sources 

The fou r sources we found with fi > 1 mas yr" 1 are all in 
the IVC94I compilation. VC94 listed four additional sources 
with fi > 1 mas yr" 1 : M87, which has a fast long-wavelength 



(18 cm) component far from the core, Cen A, which is in the 
southern sky and therefore not included in our study, and two 
others, Mrk 421 and 115 6+295, wh ere our measured values 
are well below 1 masyr" 1 dPaperllJ) . 

We note that, with years of increasingly better observations 
on more objects, the known number of sources with fast com- 
ponents has not increased. There are only 5 compact jets 
that show [a > 1 masyr" 1 at 15 GHz, within our flux den- 
sity range. These are all nearby objects and include three 
galaxies, 3C 111, 3C 120, and Cen A; one BL Lac object, 
BL Lac itself; and one quasar, 3C 273. Monthly monitoring 
at higher resolution by I J05l detected 5 sources (out of 15) that 
had [i > 1 mas yr" 1 . We found 4 of these, but we measured 
fi < 1 masyr" 1 for their fifth object, 1510-089. In addition, 
they measured fi — 1 mas yr" 1 for 0219+428 (3C 66A), but it 
has low flux density and is not in our survey. 

5.3. The Boost Exponent 

The exponent n in equation [3] controls the sharpness of the 
peak of the aspect curve. Figure [7] shows the data with three 
aspect curves for 7 = 32, with different values of n. The curves 
have been truncated at the 4% and 96% probability limits, and 
the values of L a have been adjusted so that the curves roughly 
match the right-hand side of the data. The probability was 
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calculated with equation (A15) from lVC94l as the simulation 
described in Appendix lAl uses n = 2 and has not been calcu- 
lated for other values of n. 

It is important to compare the curves with the data only in 
the region where the probability is significant. From Figure [7] 
there is no strong reason to pick one value of n over another. 
However, if n = 3, the boosting becomes so strong that strong 
distant quasars, near the peak of the distribution, have L a as 
small as the jets in weak nearby galaxies, which (we argue 
in ® are only mildly relativistic. This is unrealistic, and we 
conclude that n < 3. The value n = 2 has a theoretical ba- 
sis 02.lt and we have adopted it here. Note the large range 
in intrinsic luminosity corresponding to different values of n. 
Since n is not known with precision, the intrinsic luminosities 
have a corresponding uncertainty. 

6. THE PEAK LORENTZ FACTOR 

Beaming is a powerful relativistic effect that supplies a 
strong selection mechanism in high-frequency observations 
of AGN. Consider a sample of randomly oriented, relativis- 
tically boosted sources that have distributions in redshift, in- 
trinsic luminosity, and Lorent z factor . M ake a fl ux density- 
limited survey of this sample. I VC94I and ILM97I have shown 
that in this case the selected sources will have a maximum 
value of /3 app which closely approaches the upper limit of the 
7 distribution, even for rather small sample size. This comes 
about because the probability of selecting a source is maxi- 
mized near 9 = Q.69 c where /3 app rj O.97 (see Figure Ht); and 
there is a high probability that, in a group of sources, some 
will be at angles close to 9 C , where /3 app « 7. Hence, be- 
cause /3 ap p.max ~ 32, the upper limit of the 7-distribution is 

7app,max ~ 32. 

7. THE DISTRIBUTIONS OF 7 AND L 

We showed in ^5. II that the lack of sources to the left of the 
envelope is not a selection effect, but is intrinsic to the ob- 
jects. Since the envelope is narrow at the top, /3 app and L are 
correlated; high /? app is found only in sources that also have 
high L, but low /3 app is found in sources with all values of L. 
This translates into 7 having a similar correlation with L a , for 
the quasars. The 7 distribution will be similar to the /3 app dis- 
tribution in Figure 7, but flatter; with many points shifted up, 
but nearly all by less than a factor 2 above /3 app . The L a dis- 
tribution will remain more spread out at low 7 than at high 7, 
leading to the correlation that the highest 7 are found only in 
jets with hi gh intri nsic luminosity. This is consistent with a 
result from 1LM97L that Monte-Carlo simulations with nega- 
tive correlation between 7 and L a give a poor fit to the statis- 
tics of the flux densities from the Caltech-Jodrell Bank survey 
dTavloret al.lll996l) . 

The good fit of an aspect curve as an envelope to the data 
in Figure [6] suggests that the parameters of the curve, 7 = 32 
and L„ = 10 25 WHz -1 , reflect the peak values of 7 and L Q in 
the populati on. The distribution of 7 may be a power law, as 
suggested by LM97 and, as discussed in ^6] 7 max = 32 is close 
to the maximum value in the distribution. We now consider 
constraints on the peak value for L a . 

Figure [8] is similar to Figure [6j but with several aspect 
curves, each showing only the region 0.04 < P(9\jf) < 0.96. 
The envelope is now formed by a series of aspect curves, 
with successively lower values of 7. Most of the sources will 
have 7 rather close to /3 app , but some will have 7 substan- 
tially greater (see Figures^ and|3^). In Figure [8] these latter 
sources will not lie near the top of an aspect curve, but will 



be down from the peak. It is more likely that they will be at 
small angles (9 < 9 C ) than at large angles. 

Consider the BL Lac marked B, near the intersection of 
curves 7 = 6 and 7 = 20. It could be on either curve, but it 
is near the low-probability region of curve 7 = 20. For ev- 
ery source on curve 7 = 20 near the intersection, there should 
be several farther up the curve. Note that the 7 = 20 curve 
intersects the limit curve K to the left of the peak. There is 
little available redshift volume at the peak, but the volume in- 
creases rapidly at lower /3 app , and the lack of sources there 
means that source B is unlikely to have 7 = 20. Alternatively, 
it could be close to an extension of curve 7=15, but then it 
is again in a low-probability region. There are a number of 
sources near the peak of curve 7 = 15, and B could be a high- 
angle version of one of them. But the probability of that is 
well below 0.04, and there can be few such sources in the en- 
tire sample of 1 19. We conclude that the galaxies and BL Lacs 
on the left side of the distribution (L < 3 x 10 25 WHz" 1 ), with 
high confidence, are not off-axis versions of the powerful 
quasars (curves 7=15 and 32), nor are they high-7, low-L e 
sources (curve 7 = 20). 

Source B in Figure [H is the eponymous object BL Lac, 
2200+420. iDenn. Mutel & Marscheri (l2000h studied BL Lac 
in detail, and showed that the jet lies on a helix with axis 
9 = 9° and pitch angle 2°. If 9 = 9° ± 2° is combined with 
our value for the apparent transverse speed, /3 app = 6.6 ± 0.6, 
then 7 = 7 ± 1. This agrees with our conclusion above. 

Now consider the sources near point C in Figure [8] at 
L = 3 x 10 28 W Hz" 1 , /3 app = 9. They could have 7 « 9, but in 
that case there should be several others down the 7 = 9 curve 
to the right, where most of the probability lies. But there are 
none there. Any aspect curve with a peak farther to the right is 
unlikely to represent any of the measured points, and so curve 
7 = 9 is about as far to the right as should be considered. If the 
sources at C are on curve 7 = 9, then their intrinsic luminos- 
ity is an order of magnitude greater than that for the sources 
near the top of the distribution, the fastest quasars. To avoid 
a negative correlation between 7 and L a , some of the sources 
near point C should have 7 = 20 or more, with the appropriate 
small values of 9. However, others near the right-hand side of 
the distribution might well have 7 ~ 9 or smaller. This means 
that the distribution of L a could extend up to 10 26 W Hz -1 . 

8. QUASARS AND BL LACS WITH (3 APP < 3 

Twenty-two of the 92 quasars and 3 of the 13 powerful 
BL Lacs in Figure[8](L > 3 x 10 25 W Hz" 1 ) have /3 app < 3, and 
have low probability if 7 > 10. What are the intrinsic proper- 
ties of this group? We consider three possibilities, (i) They are 
high-7 sources seen nearly end-on, and have P(9) < 0.04. We 
expect only a few such end-on sources out of a group of 105. 
Most of the low-speed quasars cannot be explained this way. 
(ii) They are low-7 high-L sources, and have 7^3. We dis- 
cussed this above for point C in Figure [8] with /3 app = 9; now 
we are considering /3 app < 3, and the argument is stronger. 
Unless the most intrinsically luminous sources have low 7, 
this option is not viable, (iii) A more likely situation is that 
many of these low-/3 app components appear to be slow be- 
cause 7,, < 7*. 

In support of comment (iii), we note that one of the slow 
objects, 1803+784, was also observed by |J05| at 43 GHz. 
They find (3 app = 15.9 ± 1.9, whereas, at 15 GHz, we found 
/?app = ~0-6 ± 0.6. The higher resolution at 43 GHz is crucial 
in detecting fast components in sources like this, because they 
are within 1 mas of the core, at or below the resolution limit at 




FIG. 8. — The data are plotted as in Figure[6] The aspect curves are truncated at the 4% and 96% probability levels. At the peak of each curve, P(8) 0.75; 
i.e. about 3/4 of the probability of selecting a source with this value of 7 and L„ is on the right side of the curve. The cross close to /3 app = marks the source 
1803+784, and point J is the same source but with the /3 app value at 43 GHz from J05, see text. Curve K is a short section of the limit curve K from Figure[6] 



15 GHz. In Figure[8]the 15 GHz speed for 1803+784 is shown 
with a cross; and the 43 GHz /3 app value, with the 15 GHz lu- 
minosity, is shown with the letter J. It is likely that we have 
reported a component speed that is not indicative of the beam 
speed for 1803+784. 

In Figure 8, source 1803+784 is in a cluster of objects that, 
formally, have negative speed. However, they all are within 
lcr of zero, and their negativity is of little significance. A 
number of other sources have components with similar for- 
mally negative speeds, but in addition they have a component 
with a larger positive speed. In this paper we have only used 
the fastest component in each source. 

9. GALAXIES 

The points in Figure [6] appear to run smoothly from low to 
high apparent luminosity, suggesting that the different types 
of objects might be closely related. However, the smoothness 
is supplied by the BL Lac objects, which connect the galaxies 
and quasars that otherwise are widely separated in apparent 
luminosity. In addition, the galaxies all have z < 0.2, and 
nearly all the quasars have z > 0.4. The separation is at least 
partly the result of our restricted sensitivity, coupled with the 
luminosity functions. We cannot observe "galaxies" at high 



redshift because our sensitivity is too low; and we see few 
"quasars" at low redshift because their local space density is 
so low. In this section we consider whether the galaxies and 
quasars form separate classes, or, in particular, whether the 
galaxies might be high-fl coun terparts of the more luminous 
sources dUrrv & PadovanjnH!. 

Three galaxies have superluminal components, and their 
speeds place them with the lower-speed quasars, as seen in 
Figure[8] These fast galaxies, shown in Table 1, all have broad 
emission lines and are classified as Syl; they are at low red- 
shift and are highly variable at radio w avelengths. The ob- 
scurin g torus paradigm for Syl galaxies dAntonucci & Millerl 
11985b suggests that they are not at large values of 8, and this 
is confirmed by the observed values of /3 app , which show that 
6 must be less than 6» max = 2arctan/3~ 1 ~ 20° to 45°. To esti- 
mate values of 7. 9, and L Q for these galaxies, we combine the 
measured /3 app with a variability Doppler factor, D val , derived 
from t he time scale and strength of variations in flux density 
(e.gJCohen e"tai]|2004l) D v „ is given bvllol for 0415+ 379 
and 0430+052, and by iLahteenmaki & Valtaoial d 19991) for 
1845+797. We have converted the last value to the cosmol- 
ogy used in this paper, and use an intrinsic brightness temper- 
ature 7}, = 2 x 10" K. This is a characteristic lower limit for 
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TABLE 1 
Galaxies with m > 1 



IAU name Alias Type Redshift /3 app D vaI 7 9 (deg) L„ 

(1) (2) (3) (4) (5) (6) (7) (8) (9) 

0415+379 3C111 Syl 0.049 6.1±0.1 3.4 7.3 15 4.4 X 10 23 

0430+052 3C120 Syl 0.033 3.6±0.2 2.4 4.1 22 4.8 x 10 23 

1845+797 3C 390.3 Syl 0.057 2.3±0.1 0.9 3.8 42 3.1 x 10 24 



sources in their highest brightness states dHoman et al.l 2006). 
and should be more appropriate than the canonical equiparti- 
tion value, for variability measurements based on flux density 
outbursts. Note that lJ05l use a different procedure to calculate 
Z) var , and do not assume an intrinsic temperature. The D var are 
model-dependent and their reliability is difficult to assess. 

The Lorentz factors for the quasars are not estimated in this 
paper, but from Figure 3 it can be seen that many of them must 
have Lorentz factors close to their apparent speeds. Thus, 
from Table 1 and Figure 8, the Lorentz factors of the super- 
luminal galaxies are comparable with those for the slower 
quasars. Their luminosities, however, do not overlap with 
those for the quasars, indicating that they are a different pop- 
ulation. 

Only one of the seven slow galaxies (/3 app < 1) has /3 app con- 
sistent with zero {2a). The others show definite motions and 
several must be at least mildly relativistic, with /3 app > 0.3. 
The galaxies with /3 app < 0.3 cannot contain a highly relativis- 
tic jet, for that would force 9 to be unacceptably small. For 
example, if /3 app = 0.3 and 7 = 5, then 9 = 0.35° and S = 9.9. 
This gives 9/9 c = 0.03, which is extremely unlikely, for an ob- 
ject with 7 = 5 (FigureQJ). In any event, all proposed galaxy- 
quasar unifications place the galaxy at a high angle, where the 
flux density and apparent speed are reduced, but a relativis- 
tic beam cannot show /3 app = 0.3 at any angle not near 0° (or 
180°). Hence, the galaxies are neither low-angle nor high- 
angle versions of the distant quasars. However, Cygnus A 
m ay be an exceptio n , as d iscussed in the next section. 

iGiovannini et al.1 (12001 1) have concluded that most radio 
galaxies (including FRI) contain relativistic jets. They as- 
sumed that all sources have jets with intrinsic bipolar sym- 
metry, and used the measured side-to-side ratio with a corre- 
lation between lobe power and intrinsic core power to obtain 
limits on j3 and 9. Our procedure may be more robust be- 
cause each source has a measured /3 app , and we do not appeal 
to symmetry of the lobes. 

9.1. Cygnus A 

The galaxy Cygnus A (z = 0.056) merits special discus- 
sion. The radio lobes are exceptionally powerful, and their 
luminosity is comparable to that of the most powerful and 
distant radio galaxies. T he jets, however , are weak. Op- 
tical polarization studies (Og le et aT1ll997t) reveal polarized 
broad lines and sh ow that Cygnus A is a modest quasar. 
iBarfhel et al.l ([1995) used the front-to-back ratio of the jets 
of Cygnus A at 6 cm, together with /3 app , to estimate 9. We 
repeat their analysis with our value for /3 app , 0.83 ±0.12, and 
obtain 45° < 9 < 7 0° . This agrees w ith other estimates of 
the angle, including lOgle et al.1 (1 1997b who found 9 > 46°, 
and lVestergard & Barthell (1 19931) . who found 9 ~ 50°to 60°. 
The combination of /3 app and 9 gives 1.24 < 7 < 1.36 and 
0.59 < (3 < 0.68. Cygnus A is mildly relativistic. 

On the other hand, because the lobes in Cygnus A are 



so powerful, we might have expected that it would have a 
highly relativistic jet. These contradictory ideas can be recon- 
ciled with a two-component beam, consisting of a fast spine 
with a slow sheath, as s uggested by numerical simulations 
(e.g.. lAgudo et al.ll2001l) . The slow beam that we see has 
Ssiow = 1 .5 Jy. The fast beam is at a high angle to the LOS and 
not seen because it is deboosted, and it must be at least a factor 
of ten weaker than the slow beam; i.e., 5f ast < 0.15 Jy. If the 
fast beam has a Lorentz factor of about 10, then if observed at 
a small angle its flux density would be up to a few hundred Jy, 
far higher than observed in any other source. But Cygnus A is 
much closer than most superluminal sources; and if the near- 
est quasar, 3C 273, were at the distance of Cygnus A, its flux 
density would be within a factor 2 of our putative value for 
Cygnus A. If Cygnus A were at z ~ 1, and pointed near the 
LOS, it would be a normal quasar, with radio and optical lu- 
minosity somewhat below the median for quasars. It is ex- 
ceptional only because it is accidentally nearby. This model 
solves the long-standing problem of the strong lobes com- 
bined with the weak core. 

In this model the total flux density from Cygnus A varies 
more slowly with 9 than S 2 , the commonly assumed law. If 
this is correct, and if it applies generally to many other radio 
sources, then it will affect the usual discussi ons of the unifi- 
cation of radio sources by aspect. See, e.g., IChiaberge et al.1 
(2000), who invoke a two-component model in their article 
on unification. 

10. DISCUSSION 

The aspect curve in Figure [6] is a good envelope to the 
quasar data, and this suggests that the relativistic beam model 
is realistic. The slow rise and rapid fall of the curve is a direct 
consequence of Doppler boosting combined with time con- 
traction, both of which are relativistic effects. We have used 
the common assumption that the moving VLBI component is 
traveling with the beam; i.e., 7,, w 72,. But if this is incorrect, 
for example, if the beam and pattern speeds are independent, 
or often are far apart, then it is hard to see why an aspect curve 
should form an envelope. In particular, if they are indepen- 
dent, then some sources should have high j p with low 7/,, and 
this could place them to the left of the envelope in Figure [6] 
The lack of sources there has been emphasized earlier, and is 
evidence that, for the fastest components in many sources, 7^ 
and jt> are closely related. , 

In about a fourth of the BL Lacs and quasars the fastest 
component appears to be moving slowly (/3 app < 3). These 
objects are discussed in ^8] where it is shown, using probabil- 
ity arguments, that most of them cannot be high-7 sources 
with a small value of 9. However, other evidence sug- 
gests strongly that most of these objects do have a highly 
relativis tic beam. This evidence consists of rapid variabil - 
ity (e.g jAller, Aller, & Hughesl2003l;lTerasranta et ai1l2005l) ; 
high apparent brightness temperature measured both directly 
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with VLBI (Pap er iVf) and indirectly with i nterstellar scintilla- 
tions dKraus et all2003l:lLovell et alJ2003h : gamma-ray emis- 
sion (e.g.. iDondi & Ghisellinil 119951) : and observational ef- 
fects due to differ ential Dopple r boosting, such as the Laing- 
Garrington effect (lLaingi ri988). Many of these sources must 
have 7 P <C 7/,, as discussed in $8] 

11. CONCLUSIONS 

1. The aspect and origin curves provide a useful way to 
understand the relations between the intrinsic parameters of a 
relativistic beam, Lorentz factor 7 and intrinsic luminosity L , 
and the observable parameters, apparent transverse speed /3 app 
and apparent luminosity L. Limits to the intrinsic parameters, 
for a given observed source, are found on its origin curve that 
has been truncated by probability arguments. 

2. About half the sources with /3 app > 4, that are found in a 
flux-density limited survey, will have 7 within 20% of /3 app . 

3. The 2-cm VLBA survey has yielded high-quality kine- 
matic data for 119 compact radio jets. When plotted on the 
observation plane, they are bounded by an aspect curve for 
7 = 32, L = 10 25 W Hz" 1 , that forms a good envelope to the 
data at high luminosities. From this, with probability argu- 
ments, we find that the peak Lorentz factor in the sample is 
7 k, 32 and the peak intrinsic luminosity is L ~ 10 26 W Hz" 1 . 

4. There is an observed correlation between /3 app and L for 
the jets in quasars: high /3 app is found only in radio jets with 
high L. This implies a similar correlation between 7 and L„: 
high 7 must preferentially exist in jets with high L . 

5. The Doppler-boosting exponent n for a typical source 
in the survey must be less than 3, or else the highly luminous 
jets with the fastest superluminal speeds will have intrinsic 
luminosties comparable to the slow, nearby galaxies. 

6. There are too many low-speed (/3 app < 3) quasars in the 
sample, according to probability arguments. It is likely that 
some of them have pattern speeds substantially lower than 
their beam speeds. 



7. The galaxies have a distribution of Lorentz factor up 
to 7 = 7; three show superluminal motion, but most are only 
mildly relativistic. They are not off-axis versions of the pow- 
erful quasars. Cygnus A may be an exception, and we suggest 
that it might have a "spine-sheath" morphology. 

8. Our results strongly support the common relativistic 
beam model for compact extragalactic radio jets. The pattern 
and beam speeds must be approximately equal, for the fastest 
components in many sources. 
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APPENDIX 

A. MONTE-CARLO CALCULATIONS 

M. Li ster et al. (in preparatio n) are performing Monte-Carlo calculations to simulate observations being made in the MOJAVE 
survey dLister & Homanl [20051) . One of their simulations that illustrates the probability densities we need is used here. The 
calculation selects 100,000 sources with S > 1.5 Jy {S > 2 Jy in the south) from a large parent population with a distribution of 
Lorentz factors having a power law with index -1.25 and peak 7 = 32, and a power-law distribution of intrinsic luminosities with 
index -2.73 and minimum luminosity 1 x 10 24 W Hz" 1 . The model uses an evolving luminosity function based on a fit to the 
quasars in the Caltech-Jodrell Bank 6-cm survey (LM97). The calculation does not assume any correlation between the intrinsic 
quantities 7 and L a . 

Figure [9] shows 14,000 of the selected sources, in (#,7) space. The curve 7sin# = 1.0 is shown and it can be seen that the 
majority of sources have 9 < 6 C . The lines 7sin# = 0.15 (bottom) and 7sinf? = 2.0 (top) are also shown; these show the 4% and 
96% cumulative levels for the simulation. They are used in the text as practical limits. 

Lister et al. discuss the distribution functions for various intrinsic and observed quantities. Here, we only show results for 
subsamples representing slices through the full distribution at constant 7, and at constant /3 app . Figure QJ shows a slice for 
14.5 < 7 < 15.5 (N=3638); the histograms are the probability density, p(9\j « 15), and the cumulative probability, P(9\ r y « 15). 
The probability density p{9\^f) varies slowly with 7, and Figure^ shows p{sin9 < 7" 1 ), the expected fraction of sources that 
will lie within their critical angle. FiguresQJ and[2j3 show similar distributions for 14.5 < /3 app < 15.5 (N=3191). 
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